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ABSTRACT 


This paper, the first of a two-part series, presents results of an unsteady rotordynamic 
analysis of a long-bearing squeeze film damper executing circular centered orbits 
using a fluid circuit approach. A series of nodes and branches represent the 
geometry of the flow circuit. The mass and momentum conservation equations are 
solved to predict the pressure distribution in the squeeze film. The motion of the 
bearing is simulated by the variation of geometry within the flow path. The modeling 
methodology is benchmarked against published experimental long-bearing squeeze 
film damper test results. The model provides good agreement with the experimental 
damping coefficient. 

Nomenclature 

C Squeeze film damper radial clearance (difference in radii) 

C p = p(oR 2 /C 2 Pressure conversion factor 

C„ = -f/e Dimensionless direct damping coefficient 

d Damper journal diameter 
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= f/s Dimensionless direct inertia coefficient 
e Dynamic eccentricity (orbit radius) 

L Damper journal length 

m Mass flow rate 

p Local pressure within the squeeze film region 
r Radial coordinate 

R Damper journal radius 

Re = po)C 2 /|i Reynolds number 
z Axial coordinate 

e = e/C Eccentricity ratio (dimensionless orbit radius) 

0 Circumferential coordinate 

p Fluid absolute viscosity 

p Fluid density 

if Shear stress 

co Frequency of damper journal (whirl frequency) 

Introduction 

Squeeze film dampers (SFD) have been used to provide damping in high speed 
rotating machinery with rolling element bearings for years. Extensive research has 
been performed both experimentally and analytically to examine the effects of 
geometry, cavitation, Reynolds number, etc. on squeeze film damper pressures and 
forces. The experimental work of Jung et al. 1 ’ 2 examined the pressure distributions 


2 



and force coefficients for a squeeze film damper executing circular centered orbits for 
both open ended (short bearing) and partially sealed (long bearing) configurations. 
Analytical, circular-centered orbit, SFD models have been published for a variety of 
Reynolds numbers for both long bearing (axial flow neglected) and short bearing 
(circumferential flow neglected) assumptions. Analytical work by San Andres 3-6 and 
Tichy 7 ' 9 addressed the effects of fluid inertia on both short and long bearing squeeze 
film dampers, and the effect of turbulence on short bearing squeeze film dampers. 
Recently, Schallhom, et al. n , described a fluid circuit approach for modeling long 
bearing squeeze film damper rotordynamics. In this paper, a steady state modeling 
technique, utilizing a fluid network code, was described and the technique was 
benchmarked against published experimental data. That modeling technique has been 
extended to an unsteady model of the damper by Schallhom l3, ,4 . 

This paper presents the results of the unsteady modeling for a long bearing squeeze 
film damper executing circular orbits. Comparisons of the unsteady analytical 
predictions to the experimental work of Jung et al . 2 are presented. 

Fluid Flow Code 

A general purpose, one-dimensional, network flow analysis computer code was 
chosen as the platform for SFD modeling 10 . The code uses a series of nodes and 
branches to define a flow network. Nodes are positions within the network where 
fluid properties (pressure, density, etc.) are either known or calculated. Branches are 
the portions of the flow network where flow conditions (geometry, flow rate, etc.) are 
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known or calculated. The code uses a finite volume approach with a staggered grid. 
This approach is commonly used in computational fluid dynamics schemes 
(Patankar 15 , Patankar and Karki 16 ). The staggered grid approach uses overlapping 
control volumes where the conservation of mass and the conservation of momentum 
are calculated in separate control volumes. For a conventional grid where only one 
control volume is used for both scalar and vector quantities, interpolation is required 
in formulating the conservation of the scalar quantities at the edges of the control 
volume and numerical errors result. For a staggered grid, however, mass flow rates 
are available at the surface of the continuity control volume and do not have to be 
interpolated. 

Governing Equations 

The primary governing equations used in this analysis are the continuity momentum 
equations. Additionally, two auxiliary equations are necessary for closure of the 
unknown variables. These auxiliary equations consist of the energy equation and the 
equation of state. Equations 1 and 2 are the unsteady governing equations used in the 
analysis. Equation 1 is the unsteady continuity equation for the flow into a node, 

where m in is the summation of flow into a node through the branches connecting 

t° it, ^ ni out is the summation of flow out of a node through the branches connecting 

to it and m source is a mass source term (which accounts for any other mass flow into 
or out of the node). Equation 2 represents the one-dimensional conservation of 
momentum for a branch in the flow network, where Ap is the pressure difference 
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across the branch, A cross Section is the cross-sectional area, L is the flow path length in a 
segment of the damper, p is the density of the fluid, u is the fluid velocity, v is the 
normal velocity of the inner segment of the damper, A,,^, is a segmented area of the 
inner race corresponding to its normal velocity, xyis the shear stress and A shear is the 
shear area. 
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Equations 1 and 2 are solved in a simultaneous solution iterative scheme (Newton- 
Raphson scheme) in order to obtain the branch flow rates and node pressures. 


In order to obtain the mass resident in a node at a given time, an additional equation is 
necessary. The equation chosen to satisfy the unknown resident mass was the 
generalized equation of state; however, temperature is necessary for the state 
equation. In order to satisfy the requirement for temperature in the state equation, the 
energy conservation equation is used. Equation 3 is the energy equation used in the 
code, where h is the enthalpy of the fluid (from which temperature can be obtained) 
and Q Source is any additional heat source added to the fluid. Note that the term, 
P A normal v > represents the work done on the fluid by the inner segment of the damper 
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(i.e. the work done on a fluid by a moving boundary). Equation 4 is the generalized 
equation of state used in the code, where R is the gas constant of the fluid, T is the 
temperature of the fluid and Z is the compressibility of the fluid. 


^ ^in^in Qsource normal ^ ^ ’ ^out^out 



Unsteady Long Bearing Squeeze Film Damper Modeling 

The model used in the analysis is shown schematically in Figure 1 . In the analysis, 
the bearing remains wrapped; however, curvature is neglected since the clearance to 
radius ratio is small (-1/40). This model varies from prior steady state modeling 
performed by the authors 12 . The motion of the inner race of the damper is simulated 
by the variation of the geometry of the flow path, as is illustrated in Figure 2. This 
geometry variation generates a pressure variation around the circumference of the 
damper. For rotordynamics applications, the pressure distribution of the fluid circuit 
model is integrated to obtain the radial and tangential forces of Equations 5 and 6. 


F r = J/Anormal C OS0d0 

(5) 

F t = 0A normal sin 6d0 

(6) 
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In a dimensionless form, the radial and tangential forces are calculated using 
Equations 7 and 8, where Cp is the pressure coefficient defined in the nomenclature. 


f = f 2K _L cos 0 d e 
Jo c p 

f = f 2 " — sinGdB 
Jo r 


In rotordynamic models, these force components are represented by rotordynamic 
coefficients. These coefficients are analogous to the mass, damping and stiffness 
terms for a spring-mass-damper system. For an uncavitated squeeze film damper, the 
only two coefficients which occur are a radial inertia (or added mass) term, D^, and a 
circumferential damping term, C„. In order to compare with published experimental 
results 2 , the damping and inertia coefficients are normalized by the eccentricity ratio 
and pressure coefficient, equations 9 and 10. Equation 9 is the non-dimensional 
circumferential damping coefficient and Equation 10 is the non-dimensional radial 
inertia coefficient. 


C„ = -f/s 

(9) 

D rr = f/8 

(10) 
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Squeeze Film Damper Experimental Benchmark 

The experimental results of Jung, et al 2 were used to validate the modeling 
methodology. The experimental work simulated a long bearing SFD by using an O- 
ring seal on one side of the bearing, and a serrated piston ring on the other side of the 
bearing to limit axial flow. A Reynolds number of 49.0 and an eccentricity ratio of 
0.82 were chosen for comparison as the data indicates that cavitation was not present. 

Results 

Figure 3 shows the pressure profile of the benchmark squeeze film damper model 
compared with Jung’s 2 experimental results. The pressure has been non- 
dimensionalized using equation 16 for comparison with Jung’s published data. 


P _ P _ P 

Preference Cp Re P R 2 © 2 


( 11 ) 


The predicted pressure profile follows a similar trend as the experimental data. The 
predicted pressure profile deviates from the experimental results in the low pressure 
region; however, the model is able to predict the phase shift observed in the 
experimental benchmark, unlike an earlier published steady state model 12 . The 
rotordynamic coefficients for the benchmark model are presented in Table 1, along 
with Jung’s experimentally derived coefficients. As Table 1 indicates, the damping 
coefficient predicted using the unsteady model lies just inside of the uncertainty band 
of Jung’s experimental data. It should be noted that the uncertainty analysis derived 
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in reference 13 is very conservative (in favor of Jung, et al .); therefore, the unsteady 
model prediction is most likely to be well within the actual uncertainty band if all 
error sources were accounted for. The unsteady model dramatically over predicts the 
inertia coefficient compared to Jung’s experimental results (unlike the steady state 
model which substantially under predicts the inertia coefficient). Figure 4 
demonstrates the periodic nature of the model. The pressure profile shown in Figure 
4 follows the variation in the clearance during the orbit (Figure 2); therefore verifying 
that for circular centered orbits, any snapshot in time should produce the same 
rotordynamic coefficients. 

Parametric Studies 

A set of parametric studies was performed on the unsteady model to observe the 
behavior of the model under various conditions for circular centered orbits. The first 
parametric study examined the variation of the Reynolds number for a fixed orbit 
size. The second parametric study examined the variation of orbit size for various 
Reynolds numbers. 

Parametric Study on Reynolds Number for Fixed Orbit Size 

In order to compare the unsteady model with the results of the benchmark 
experimental work 2 , a parametric study on the prediction of rotordynamic coefficients 
over a range of Reynolds numbers was performed. This study was carried out for the 
following Reynolds numbers: 15, 25, 35, 49 and 75. Figure 5 shows the 
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dimensionless pressure profiles for Reynolds numbers 15, 25, 35, and 75, 
respectively. 

Figures 6 and 7 show the coefficients predicted by this parametric study compared to 
Jung’s experimental data and his prediction. In Figure 6, the unsteady model over- 
predicts the damping coefficient by ~ 10-15%, when compared with the test data 
(considered excellent agreement in rotordynamics). It is noted, however, that the 
experimental data has significant scatter and the precision uncertainty has not been 
included in the uncertainty estimate. The model prediction of the damping coefficient 
is an improvement over Jung’s analytical prediction. As Figure 7 shows, the unsteady 
model over-predicts the inertia coefficient compared with the experimental data. The 
reason for this over prediction is the difference between the predicted and 
experimental pressure at the position of lowest pressure (see Figure 3 for 0 < 0 < 1). 

Parametric Study on Orbit Size for Various Reynolds Numbers 

A parametric study was performed to examine the effect of orbit size (orbit radius) on 
the prediction of rotordynamic coefficients. Three orbit eccentricity ratios (e) were 
chosen for the study: s = 0.25, 0.50, and 0.82, respectively. Recall that 8 = 0.82 is the 
eccentricity ratio of the benchmark case. Figures 8 and 9 are the dimensionless 
pressure profiles for eccentricities of 8 = 0.25 and 8 = 0.50 at three Reynolds 
numbers: 15, 49, and 75, respectively. Figure 10 is an overlay of the pressure profile 
for the three orbit sizes at a Reynolds number of 49. As orbit size decreases, the 
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circumferential pressure variation reduces and becomes less skewed (the curvature on 
the low pressure end resembles that of the high pressure end). As a result, the 
damping coefficient and force decrease with decreasing orbit size. 

The predicted dimensionless rotordynamic coefficients are shown in Figures 1 1 and 
12. Figure 11 illustrates that the dimensionless damping coefficient increases with 
increasing orbit size and Reynolds number. In Figure 12, the dimensionless inertia 
coefficient decreases with increasing orbit, the opposite of damping. However, in the 
definitions of these coefficients, the eccentricity ratio appears in the denominator (see 
Equations 9 and 10). It is this division by £ which causes the reversal in order 
between Figures 11 and 12; therefore two additional figures (Figures 13 and 14) are 
provided to illustrate the non-dimensional tangential and radial forces (which 
correspond to the damping and inertia forces, respectively, see Equations 7 and 8). 

Figure 13 further illustrates that the tangential (i.e. damping) force increases with 
increasing orbit size. This result is expected based upon rotordynamic theory 6 . In 
Figure 14, however, the effect of orbit size on the radial (i.e. inertia) force is 
inconclusive, since the orbit size does not appear to correlate with the radial force. 
This lack of correlation is most likely due to the over-prediction of the low pressure 
region, as was seen in the comparison between the unsteady GFSSP prediction and 
experimental data for the benchmark case (Re = 49, £ d = 0.82). 
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Conclusions 


An unsteady long bearing squeeze film damper (SFD) rotordynamic analysis has been 
developed and successfully implemented. The approach uses a general one 
dimensional fluid flow code in which shear stress is calculated using friction factors. 
The motion of the rotor is simulated by variation of the geometry of the flow circuit 
and accounting for the affect of a cross inertial term in the momentum equation due to 
this geometry change. Results show that the unsteady modeling technique works over 
a wide variety of conditions. Reynolds numbers ranging from 15 to 75, and 
eccentricity ratios ranging from 0.25 to 0.82 were examined. Results of a benchmark 
model provided good agreement with experimental data (Jung, et ai, 1991b) for the 
prediction of the damping coefficient. The unsteady prediction of the damping 
coefficient varied less than 14% from experimental results, when experimental 
uncertainty is neglected. The unsteady damping coefficient results are within the 
error band when uncertainty is accounted for. The inertia coefficient is over-predicted 
by a factor of ~2 compared to the benchmark case. 

A parametric study using the benchmark geometry for a variety of Reynolds numbers 
has been performed. This study shows that the damping coefficient prediction falls 
within the uncertainty band of Jung’s experimental data, an improvement over Jung’s 
own analytical predictions. 

A second unsteady model parametric study was performed to observe the damper 
performance at various eccentricity ratios for a range of Reynolds numbers. The 
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variation of the eccentricity ratio has a large effect on the predictions of the damping 
(tangential) force component. As the eccentricity ratio increases, the damping force 
increases. However, the effect on the inertia (radial) force component is inconclusive. 
Further experimental investigation is necessary to determine the radial force variation 
with eccentricity ratio and confirm the predictions of the damping force. 

Although the predicted radial and tangential forces may be comparable in magnitude, 
a typical rotor system will have many other sources of radial force (probably of 
greater magnitude than that of the damper); however, the system may not have any 
other source of tangential force approaching the magnitude provided by the damper. 

In conclusion, it must be noted that the principle purpose of a squeeze film damper is 
to provide a damping (tangential) force to a rotordynamic system, not an inertial 
(radial) force. An unsteady analysis which can accurately predict the damping force 
over a variety of operating conditions is a valuable tool for rotordynamic analysis. 
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Table 1 : Dimensionless Benchmark Rotordynamic Coefficient Comparison 
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Figure 4: GFSSP Predicted Pressure Profile vs. Node Over One Period for Benchmark 

Case 











Figure 6: Unsteady GFSSP Model Dimensionless Damping Coefficient Predictions 




Figure 7: Unsteady GFSSP Model Dimensionless Inertia Coefficient Predictions 





Figure 8: Dimensionless SFD Pressure Distribution vs. Angle for s = 0.25 




Figure 9: Dimensionless SFD Pressure Distribution vs. Angle for e = 0.50 




Figure 10: Dimensionless SFD Pressure Distribution Comparison vs. Angle for Re = 49 





Figure 1 1 : Unsteady GFSSP Model Dimensionless Damping Coefficient - Parametric 

Study Results 
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Figure 12: Unsteady SFD Model Dimensionless Inertia Coefficient - Parametric Stud 

Results 




Figure 13:Unsteady SFD Model Dimensionless Tangential Force - Parametric Stud 
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Figure 14: Unsteady SFD Model Dimensionless Radial Force - Parametric Study Results 






